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ABSTRACT 

In this paper we report on the two-year-long Fermi-LAT observation of the 
pecuhar blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state 
from the start of science operations of the Fermi Gamma-ray Space Telescope in 
2008 August until 2009 September, and then became more active, with gradu- 
ally increasing flux and some moderately-bright flares. In 2010 April and June, 
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4C +21.35 underwent a very strong GeV outburst composed of several major 
flares characterized by rise and decay timescales of the order of a day. During 
the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law form 
with spectral breaks observed near 1 — 3 GeV photon energies. We demonstrate 
that, at least during the major flares, the jet in 4C +21.35 carried a total ki- 
netic luminosity comparable to the total accretion power available to feed the 
outflow. We also discuss the origin of the break observed in the flaring spectra 
of 4C +21.35. We show that, in principle, a model involving annihilation of the 
GeV photons on the He II Lyman recombination continuum and line emission of 
"broad line region" clouds may account for such. However, we also discuss the 
additional constraint provided by the detection of 4C +21.35 at 0.07 — 0.4 TeV 
energies by the MAGIC telescope, which coincided with one of the GeV flares of 
the source. We argue that there are reasons to believe that the < TeV emission of 
4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not likely 
to be produced inside the broad line region zone of highest ionization (~ 10^'' cm 
from the nucleus), but instead originates further away from the active center, 
namely around the characteristic scale of the hot dusty torus surrounding the 
4C +21.35 nucleus (~ 10^^ cm). 

Subject headings: galaxies: active — galaxies: jets — quasars: general — quasars: 
individual(4C +21.35) — gamma rays: galaxies — radiation mechanisms: non- 
thermal 



1. Introduction 



Radio-loud active galactic nuclei (AGN) with relativistic jets oriented at small viewing 
angles — called 'blazars' — are known for their particularly intense and dramatically vari- 
able 7-ray emission within the GeV photon energy range. This emission, strongly Doppler- 
boosted in the observer frame, dominates the total observed radiative output of a typical 
luminous blazar, e.g., flat spectrum radio quasar (FSRQ), by orders of magnitude, as estab- 
lished by the pre vious studies with the EGRET instrument onboard the Compton Gamma- 
Ray Observatory (IHartman et al.lll999[ ). There is a need for careful investigation of the 7-ray 
properties of powerful blazar sources, enabled recently by the excellent perform ance of the 

1, even 



Large Area Telescope (LAT) onboard the Fermi satellite (lAtwood et al.l 12009). Sti 



though hundreds of FSRQs have already been detected by LAT (lAbdo et al. 



spectra of only several of the brightest ones can be studied in detail ( lAbdo et al. 



2010 



b), GeV 



2010eh . 



Moreover, until now only a few luminous (quasar-hosted) blazars have been detected at 
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higher, TeV photon energies by the modern ground-based Cherenkov telescopes. 



4C +21.35 (also k nown as PKS 1222+216, z = 0.432: lOsterbrock fc Poggd Il987l) is a 
7-ray-emitting FSRQ ( iHartman et al.lll999l : lAbdo et al.ll2010a[ ). Its peculiar large-scale (~ 
100 kpc) radio structure is reminiscent of a 'Wide- Angle- Tailed' morphology characteristic 
of radio galaxies of intermediate po wer located at the centers of luminous clusters, but is 
quite unexpected for quasar sources (ISaikia et al.lll993l ). At smaller (mas) scales apparent 
superluminal velocities have been detected for a few sub- components of t he 4C +21.35 jet, 
with Papp > 10 for the Hubble constant Hq = /i lOOkms"^ Mpc~^ (IHooimeyer et al. 



19921 : iJorstad et al.l l200ll : iHoman et al.l 120011 ) . The observed superluminal blobs exhibit 
non-radial mo t ions, involving changes in the position angles and perpendicular acceleration 
(IHoman et al.l l200ll ) . The unresolved self-absorbed radio cor e is characterized b y a fiat 
spectrum and a brightness temperature of T;, ~ 5 x 10^1 K jKovalev et al.ll2005h . which 
again indicates a small jet inclination and relativistic beaming effects. Quite surprisingly for 
a blazar source, the radio core of 4C +21.35 is less luminous than its large-scale structure: 
this object is formally a 'lobe-dominated quasar', sinc e the ratio of the core-to-extended 



radio fluxes at GHz frequencies is of the order of unity (IKharb fc Shastril |2004J : IWang et al. 
200l). 



4C +21.35 was d etected by EGRET (as the most likely counterpart of 3EG J1224+2118; 
Hartman et al.l 119991 ), and subsequently subje cted to several 7-ray studies and multiwave- 



length modeling (e.g.. iNandikotkur et al.ll2007l ). During the first three months of the Fermi 



missi on, 4C +21.35 was n ot bright enough to be included in the Fermi-LAT Bright Source 
List (lAbdo et al.M2009bl ). It was detected, however, in the First LAT Catalog (IFGL; 
Abdo et al.ll2010al ). where it was associated with IF GL J1224. 7+2121. A 7-ray fiare from this 
source was noted in 2009 April (ILongo et al.ll2009t ). It was followed in 2009 December by an 



even larger fiare, seen by b oth the AGIL E Gamma-ray Imaging Detector (jVerrecchia et al. 



20091) and the Fermi LAT (ICiprinil 120091 1. On 2010 April 24 , Fermi LAT detected a par- 
ticularly stron g GeV outburst from the object (jDonatdl2010l ). The analysis performed by 
Neronov et al.l (l2010al Jbl) indicated that during the fiare the 7-ray emission of 4C +21.35 ex- 
tended up to observed photon energies greater than 100 GeV, i.e., up to the very-high-energy 
(VH E) band. At that tina e, accompanying brightening was reported at near-infrared frequen- 
cies dCarrasco eraPboioh. A sec o nd huge GeV outburst was recorded by Fermi and A GILE 
in 2010 June Jiafrate et allboiol : Istriani et aPboioh . At the same time (2010 June 17), a 
prominent excess around ~ 0.07 — 0.4 TeV photon energies w as noted by the MAGIC at 



mospheric Cherenko v telescope from the position of 4C+21.35 (lMariottill2010l : lAleksic et al. 



submitted to ApJLl ). establishing this source as the third FSRQ firmly detecte d in the VHE 



range by ground-based i nstruments besides 3C 279 and PKS 1510—089 (see lAlbert et al. 



20081 : IWagner et al.ll2010l . respectively). Followup observations by ground-based optical tele- 
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scopes ( iDominici et al.ll2010l: iNesci &: Montagnill2010[ ) and in the X-ray/ultraviolet domain 
by the Swift satelhte (IVerrecchia et aL illustrate the broad multiwavelength interest 

in the flaring activity of this blazar. 

In this paper, we study the temporal and spectral evolution of 4C +21.35 in 7 rays, as 
observed by LAT in the first two years of the Fermi operation. The paper is organized as 
follows. In § 2, we describe Fermi-L AT observations and data reduction, as well as the main 
observational findings. A general discussion is presented in § 3, while the main conclusions 
are summarized in §4. A standard ACDM cosmology with = 0.73, Qu = 0.27, and 
Ho = 71kms~^ Mpc~^ is assumed throughout the paper. 



2. Observations and Results 

The LAT onboard Fermi is an electron-positron pair-conversion telescope sensitive to 7 
rays with energies in the range from 20 MeV to more than 300 GeV. The instrument is made 
of a high-resolution silicon micro-strip tracker, a Csl hodoscopic electromagnetic calorimeter, 
and a segmented pla stic scintillator detector to identify the background of charged particles 



(lAtwood et al.ll2009l ). Compared with its predecessor EGRET, the LAT has a larger field 
of view (~ 2.4 sr), a larger effective area (~ 8000 cm^ for > 1 GeV on-axis photons), and 
an improved angular resolution {6qs ~ 0.6° at 1 GeV for events in the front section of the 
tracker). The entire sky is observed every 3 hours in a survey mode. 

The LAT data presented here were collected from 2008 August 4 to 2010 August 4. Only 
events with energies greater than 200 MeV were selected to minimize the systematic uncer- 
tainties. We used only "Diffuse" class events, which have the highest probability of being 
photons. We excluded the photons arriving from zenith angles > 105° to limit contamination 
from Earth limb 7 rays. The collected LAT data were analyzed using the standard Science 
Tools package (v9rl6pl). The Instrument Response Functions (IRF) P6_V3_DIFFUSE were 
used. To produce the hght curves and for the spectral analysis the standard tool gtlike was 
used. The photons were extracted from a region of interest (ROI) centered on the source and 
having a radius of 10°. The gtlike model (source model) includes the 4C -1-2 1.35 point source 



and all the point sources from the first LAT catalog ( lAbdo et al.l l2010af ) that fall within 
15° of the source. The model also includes a background component of the Galactic diffuse 
emission (gll_ieiii_v02 . f it) and an isotropic component (isotropic_iem_v02.txt). The 
isotropic component includes both the contribution from the extragalactic diffuse emission 
and from the residual charged-particle backgrounds. The estimated systematic uncertainty 
of the flux is 10% at 100 MeV, 5% at 500 MeV, and 20% at 10 GeV. 
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Figure 1 shows the weekly (7-day-bin) 7-ray hght curve of 4C +21.35 (at photon energies 
E > 100 MeV, extrapolating from the measurements above 200 MeV using a best-fit power 
law) since 2008 August 4 (MJD 54682). The flux F>iooMeV was nearly constant until a 
moderate flux enhancement was seen around MJD 55100 (2009 September 26). We define this 
interval as a 'quiescent state' (see also Table 1). After the enhancement, the flux gradually 
increased with some minor flares superimposed until the source entered a particularly active 
phase around MJD 55310 (2010 April 24). We define the period MJD 55100-55310 as 
an 'intermediate state', and the following one as an 'active state'. Figure 2 (a) shows the 
daily light curve (again above 100 MeV) starting from MJD 55277 (2010 March 22). As 
shown, two prominent outbursts were detected during the active state around MJD 55315 
(2010 April 29) and 55365 (2010 June 18). In the figure, black arrows indicate the epochs 
when two VHE photo ns with energies of 149 and 129 GeV were detected by Fermi-LAT 



(iNeronov et al.l l2010al ). A red arrow depicts the date MJE ) 5536 4.9 (2010 June 17) when 
MAGIC detected VHE emission from 4C +21.35 dMariottil lioiol ) : this coincides with the 
second huge outburst of the source at GeV energie^. Interestingly, during the first major 
flare the GeV flux shows "double-horn" structure; namely, a sharp rise by a factor of ~ 4 
within 1 day followed by a gradual 3-day decrease, then a gradual 3-day rise followed by a 
rapid decrease by a factor of ~ 6 within 1 day. The second flare shows the same pattern 
as the beginning of the first flare, a fast rise by a factor of ~ 3 within 1 day was followed 
by a gradual 4-day decrease. Figure 2 (b) exhibits the time variation of the photon index 
F derived by fitting the data using a single power-law function. In Figures 2 (c) and 2 (d), 
we plot the daily changes in the photon indices as a function of the 7-ray fluxes during the 
two flaring states for the two major outbursts. As shown, during the second flare — the 
one characterized by the asymmetric profile with the decay timescale longer than the rise 
timescale — a clock- wise evolution on the F — -F>iooMeV plane can be noted, while during the 
first flaring epoch (composed of two sub- flares), a hint of 'harder-when-brighter' behavior 
can be found. 

We performed a detailed spectral analysis for the accumulated dataset, investigating the 
spectral energy distribution (SED) of 4C +21.35 within the LAT energy range during each 
activity state of the source defined in Figure 1 and 2 (a). The resulting spectra are shown in 
Figure 3. As shown, the spectra are in general harder during the active states when compared 
to the quiescent state. The quiescent SED is well described by a single power-law function 



^We note that, in principle, the detection of the > 100 GeV emission from the cosmologicaUy distant 
objects (especially as distant as the source discussed here with z = 0.432) may provide interesting or even 
crucial constraints on the uncertain lev el of the Extragal actic Background Light in the infrared-to-UV 
photon energy range (see in this context lAbdo et ahluOlOcT ). This issue is however beyond the scope of our 
paper. 
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with photon index T = 2.57 ± 0.07. This value is in agreement with that given in IFGL 
for this object, and is also a typical v alue for the GeV photon indices of FSRQs detected 



by the Fermi LAT f lAbdo et al.ll2010bl ). During the intermediate state, a single power-law 
model is significantly rejected compared to a broken power-law model (the difference of the 
logarithm of the likelihood fit^ —2AL = 13.2, which can be related to the distribution 
using Wilks' theorem), consistent with a presence of a spectral break around i?br — 3 GeV. 
Similarly, during the first flare a broken power-law model is favored over a simple power-law 
(— 2AL = 35.8). The same is true for the second flare, albeit at a lower significance level 
(-2AL = 12.6). 

Table 2 summarizes the broken power-law fits (low- and high-energy photon indices Fle 
and Fhe, respectively) for the discussed epochs. It is important to investigate, however, if 
the curvature of the flaring GeV spectrum in 4C +21.35 is indeed as sharp as implied by 
the broken power-law fit, or if it is instead much smoother, being better characterized by, 
e.g., the log-parabola. The corresponding 'alternative' log-parabola fits are therefore given 
in Table 3. It can be seen that the log parabola is also preferred over the simple power 
law, though the values of AL are not as large as in Table 2. To investigate whether the 
broken power law is a better fit to the data than the log parabola we use a criterion 
based on the fluxes in 12 energy bands (see Figure 3). By combining the whole data set 
accumulated during both the intermediate and active periods (MJD 55086-55371), we found 
that the broken power-law model gives x^/d-o.f = 10.9/8, fully compatible with the data, 
while the log-parabola gives /d.o.f = 27.7/9, corresponding to a chance probability of 
~ 0.001 (equivalent to ~ 3.3o"), and is hence disfavored by the data. We therefore conclude 
that the broken-power-law model is a better representation of the source spectrum during 
the active periods. 

The three different SEDs of 4C +21.35 corresponding to different periods of the active 
state of the source are compared in Figure 4 with the flaring spectrum of another famous 
FSRQ detected by the Fermi LAT, 3C 454.3 ('th e brightest flare spectru m' of 3C 454.3 cor- 



responding to the epoch MJD 55166-55173; see lAckermann et al.ll2010l ). This comparison 
indicates a similar break energy located around 1 — 3 GeV, but a variety of low- and high- 
energy spectral slopes. Note, finally, that during the flares the 7-ray continuum of 4C +21.35 
extends up to GeV photon energies with rather flat photon indices F ~ 2, and this is a quite 
unexpected behavior for FSRQs (see the discussion below). 



^— 2AL = — 2 log(LO/Ll), wher e LP and LI are th e maximum likehhood estimated for the null and 
alternative hypothesis, respectively (jMattox et al.lll996l ). 
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Discussion 



3.1. Source Energetics 

With the assumed cosmology, the redshift of 4C +21.35 corresponds to a luminosity dis- 
tance of di^ = 2.4 Gpc. At this distance, the 7-ray spectrum of the source extending during 
the flares with flat photon indices F ~ 2 up to the observed photon energies of 1 — 2 GeV, 
where it reaches monochromatic energy fluxes of the order E Fe — 10~^ ergcm"^ s~^, implies 
isotropic GeV luminosity ATid^ x [E Fe\ ~ 7 x 10^''ergs~^. Anticipating a bolometric correc- 
tion of the order of a few, the corresponding ('flaring') 7-ray luminosity of 4C +21.35 would 
therefore be ~ lO^^ergs"^. The analogous value for the quiescent state should then be 
one to two orders of magnitude lower (see Figure 3). Also, since the 7-ray emission of FSRQs 
— generally thought to result from inverse-Comptonization of ambie nt photon fields within 
relatiyistic jets at sub-p c/pc distances from the central engine (e.g.. iGhisellini fc Tavecchio 
2OO9I : ISikora et al.ll2009l . and references therein) — is beamed in the observer rest frame, the 
total power emitted during the major flares is 



L ~ 



^ jet -^7 — 



10 



46 



Ljct 

10 



erg s 



(1) 



where Fjet is the jet bulk Lorentz factor (ISikora et al.l 119971). Note that in the broad-band 
modeling of blazar sources detected by EGRET. ICelotti &: Ghisellinil ( 120081 ) claimed Fjet — 15 
for 4C +21.35. It is interesting to compare the emerging value of Lem with the total power 
radiated by the central engine, i.e., by the supermassive black hole/accretion disk system in 
this object. 



For 4C +21.35, iFan et al.l ( l2006l ) give the energy flux of the H/? emission line Fh/3 — 
3.1 X 10~^^ erg cm~^ s~^. This corresponds to an isotropic luminosity of Lh^ ~ 2 x 10'^^ erg s~^, 
and hence to the total luminosity of the broad-line region (BLR) of Lrlr — 25.3 x Lh/j — 



5 X 10^^ ergs ^ (jWang et al.l l2004l : iFan et al.l l2006l ). This value is in a good agree ment 



with (a factor of 2 lower t han) the corresponding one claimed by ICelotti et al.l (119971 ) and 
Celotti &: Ghisellinil (120081 ) for the source. The luminosity of the accretion disk illuminating 
and ionizing the BLR clouds in 4C +21.35 can thus be evaluated as, roughly. 



disk 



BLR-^BLR 



0.5 X 10^^ 



0.1 



ergs 



(2) 



where ^blr is the fraction of the disk radiation reprocessed in the BLR (which is expected 
to be at the level of 10%). This implies that, during the 7-ray flares, the radiative output 
of the nuclear jet in 4C +21.35 is of the same order of magnitude as that of the central 
engine, Lcm/-^disk ~ (rjct/10)~^ (^blr/0.1). With the conservative estimates for radiative 
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efficiencies of an AGN jet and of a standard accretion disk, the energetics of 4C +21.35 
therefore becomes extreme, as argued below. 



First, let us note that IWang et al.l ( 120041 ) estimated the mass of the supermassive black 



hole in 4C +21.35 as A^bh — 1-5 x 10 Mq. This gives the Eddington luminosity of L^dd = 
47r GA^BH ^pc/dT — 2 X lO^^ergs"^ and the accretion-related luminosity in Eddington units 
Ldisk/-^Edd — 0.3 (^blr/0.1). The implied high accretion rate indicates the accretion disk 
in 4C +21.35 to be of the standard "optically-thick, geometrically-thin" structure, radiating 
with the efficiency of ?7disk = -^disk/-^acc — 0.1, where Lace is the total power of the accreting 
plasma. The analogous radiative efficiency of a nuclear AGN jet, which could be defined 
as 77jet = Lem/ Ljet foY the total jet kinetic power Ljet, is a poorly known parameter, though 
there are arguments in favor of r/jet ^ 0.1 at most, since AGN jets in luminous sources are 
expected to transport (and to depo sit) the bulk of their kin etic energies far away from their 



bases (on kpc-Mpc scales; see, e.g.. lSikora fc Madejskill2000l and references therein). Hence 



the efficiency of the jet production in the discussed object appears to be very high, namely 

' ' ' ~ 1 , (3) 



Lace V 10 y V 0.1 / V 0.1 / VO.l 

indicating that, at least during the flares, the jet in the analyzed source carries away a total 
power comparable to the total accretion luminosity available to feed the outflow. As a result, 
the two major flares detected by Fermi-LAT should be considered as being saturated at the 
maximum level. In other words, we should not expect to detect any flares from 4C +21.35 
more luminous than -F>iooMeV ~ lO^^phcm"^ s~^. 



It is interesting to note in this context that lAbdo et al.l (l2010dl ) — who modeled the 
broad-band spectral energy distributions of several FSRQs detected by the Fermi LAT in 
a quiescent state and compared the emerging model parameters with the analogous ones 
obtained during the flaring epochs of the same objects — concluded that "the difference be- 
tween the low- and high-activity states in luminous blazars is due to the different total kinetic 
power of the jet". This, along with our finding regarding the energetics of 4C +21.35 as pre- 
sented above, could therefore suggest an interesting scenario for the nature of the variability 
in luminous blazar sources. Namely, the two findings, when considered together, would be 
consistent with the situation in which during the quiescent states the total kinetic power of 
quasar jets constitutes only a small fraction of the accreting power (Ljet/Laec ~ 0.01 — 0.1), 
while during the major outbursts — triggered most likely by some sort of instabilities oper- 
ating within the accretion disk — the jets carry away almost all of the available accretion- 
related luminosity (Ljet ~ -^aee)- Since for the quasar-hosted blazars the fraction of time in 
which a source is in a flaring state (i.e., when the observed flux is ab out one order of magn i- 



tude or more higher than the average flux of a source) is about 1% (ITavecchio et al.ll2010l ) 



- 10 - 



such epochs characterized by the extremely efficient jet production are not representative of 
the whole kinetic output of the central engine, however. 

Even though 4C +21.35 may be considered as a peculiar blazar for a number of reasons 
(comparable core and extended radio fluxes, or a relatively small A^bh, for example), we 
believe that the above finding/speculation should apply to most FSRQs. That is because 
in terms of the standard radio-loudness parameter, this object is quite representative for 
its class. In particular, for the S-band luminosity Lb — iq^-'^\^b\ +35.6 gj^gg-i 10^^ e rgs~^ 
following from the absolute magnitude Mb — —23.92 as given in I Wang et al.l (12004), an d 



ext 



the 5 GHz radio luminosity L|j 



10 ergs as given in 



Kharb fc Shastril (120041 ) . 

10^ 



one can find the radio-loudness parameter of 4C +21.35 sour ce TZ = 10^ x [Lr/Lb 
being typical for jetted quasars in general (ISikora et al.l 120071 ). 



3.2. 7-ray Emission of 4C +21.35 



As described in the previous section (§ 2), the flaring spectra of 4C +21.35 exhibit clear 
breaks at observed photon energies 1 — 3 GeV. Such breaks, which seem characteristic for 
bright FSRQs in general, were studied most extensively in the case of the anal ogous blazar 



3C 454 .3. In particular, during the early Fermi-LAT observations of 3C 454.3, lAbdo et al. 



(l2009al ) reported on the GeV flare characterized by the rise and decay timescales of > 3 days, 
which clearly displayed a broken power-law form of the emission continuum with the low- 
and high-energy phot on indices Ft.e — 2 - 3 and Thf, — 3.5, respectively, and the break photon 
energy Ey^j- ~ 2 GeV. lAckermann et al.l (l2010l ) analyzed the subsequent major outburst of 
3C 454.3, demonstrating that despite significant flux changes in 7-rays (involving multiple 
flares with different profiles and flux doubling timescales possibly as short as hours), the 
photon break energy remained stable. Those findings are therefore compatible with the ones 
presented in our paper. 



Abdo et al.l (l2009al ) argued that the break in the GeV spectrum of 3C 454.3 cannot re- 



sult from simple cooling effects, not only because the spectral change is larger than the one 
following from the simplest h omogeneous ern ission models (AF = Fhe — Fle — 1-2 > 0.5; 
although see in this context iReynoldd 120091 ). but also because the location of the break 
— similar for different flares and different activity levels of the source — is not consistent 



with the expected one (around the observed MeV photon energies; see ISikora et al.ll2009l ). 
The latter argument holds also for 4C +21.35, even though in this case the relatively mild 
spectral chang e observed, AF ~ .5 (see Table 2), is seemingly suggestive of the standard 



cooling effects. I Abdo et al.l (l2009al ) showed also that the broken power-law form of the GeV 
continuum in 3C 454.3 cannot result from the absorption of the emitted 7-rays via photon- 
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photon annihilation on the lower-frequency jet emission, or on the extragalactic background 
light, since in both cases the corresponding optical depths for 1 — 3 GeV photons are signif- 
icantly less than unity. Finally, it was argued that the pair-production processes involving 
the direct X-ray emission of the accretion disk and disk corona cannot play a role in this 
context as well, because the blazar emission zone would then have to be located too close 
to the black h ole (a few Schwarzch ild radii). Again, the same is true for 4C-t-21.35. For all 
those reasons, lAbdo et al.l (l2009al ) advocated for the observed break in the 7-ray spectrum 
of 3C 454.3 to reflect a break in the underlying energy distribution of the radiating electrons, 
which is "intrinsic" to the particle acceleration processes taking place within the emission 
region. 



More recently, iFinke fc Dermerl ( 120101 ) proposed an alternative explanation for the origin 
of the break in the spectrum of 3C 454.3, discussing a particular model for the broad-band 
emission of the source. In the model, a break around the observed GeV photon energies is 
due to a transition from the dominant inverse-Compton scattering (in the Thomson regime) 
of the direct emission from the accretion disk, to the dominant inverse-Compton scattering 
(taking place in the Klein-Nishina regime) of the disk emission re-processed in the BLR. The 
change in the dominant seed photon population is crucial, since the transition between the 
Thomson and the Klein-Nishina regimes in the inverse-Compton scattering process involving 
a single target photo n field cannot result in a relatively sharp break in the emission spectrum 
(see the discussion in lAckermann et al.ll2010l . and refe rences therein). Howeve r, we note that 
in evaluating the 7-ray output of the 3C 454.3 jet, iFinke fc Dermerl (120101 ) did not take 
into account any line emission of BLR clouds ionized by the accretion disk, but considered 
only a disk continuum Thomson-scattered by the BLR clouds as a soft photon field for the 
inverse-Compton radiation; in reality, the line emission may dominate over the Thomson- 
scattered continuum by at least one order of magnitude (e.g. iTavecchio fc Ghisellinil |2008| . 
and references therein), and this can possibly affect the model results. 



This line emission, on the other hand, was analyzed carefully by iPoutanen fc Stern 

: n their attempts to explain the GeV breaks of bright Fermi blazars. In particular, 
Poutanen fc SternI ( 120101 ) argued that the pair production on the He II Lyman recombination 
continuum and line emission of BLR clouds may be responsible for broken power-law 7-ray 
spectra of FSRQs, with spectral breaks formed around photon energies of ~ 5 GeV. The 
model, if correct, would imply the blazar emission zone (at least in the case of quasar 
sources) to be located "inside the region of highest ionization of the BLR", i.e. at < 0.1 pc 
distances from the central engine. One should mention at this point that a possible challenge 
to the proposed scenario could be a lack of any observational signatures for the emission of 
secondary pairs created via the photon-photon annihilati on in the spectra of FSRQs (e.g., 
in the X-ray domain; see in this context the discussion in iGhisellini &: Tavecchidl2009l ). but 
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this issue would have to be investigated in a more quantitative way in the particular context 
discussed here before drawing any definitive statements. In addition, in the model one 
should also expect the Klein- Nishina effects to modify the energy distribution of the highest- 
energy electrons cooling radiatively via inverse- Compton upscattering of the UV photons 
from the BLR, an d thus producing the obs ervable signatures in the synchrotron continua of 
l uminous blazars (Moderski fc Sikorall2005l ). Before now, no such signatures have been found 
( ISikora et al.ll2009l ). although the frequency regime of interest (UV/soft X-ray) is particularly 
challenging observationally. Keeping the above cavea ts in mind, below we discuss in more 
detail whether the model by lPoutanen &: SternI (120101 ) can be applied to 4C +21.35. 



First, let us therefore estimate the luminosity of the He II Lyman complex in 4C +21.35. 
For a wide range of the values characterizing the ionization parameter and column den- 
sity in the closest vicinities of black hole/accretion disk sys tems of quasar sources, one 



expects I/HeTT T.y to be of the order of 10% of the Ly a emission ( jTavecchio fc Ghisellinill2008 



Poutanen fc SternI l2010l ). Hence, for the Lya luminosi ty of the BLR in 4C +21.35, namely 



4.5 X L 



H/3 



lO^'^ergs ^ 



,see 



Wang et al.1 120041 ) . we get Lneii Ly — 10"^^ ergs ^ Sec- 



ond, the characteristic scale for the BLR in the object, -Rblr, should be evaluated. This can 
be done by applying the stan dard relation between the disk luminosity Ldisk and -Rblr (e.g., 
Ghisellini fc Tavecchioll2008l . and references therein), to obtain 



R 



BLR 



0.1 



La; 



disk 



1046 erg g-i 



1/2 



pc 2 X 10 



17 



BLR 



-1/2 



cm 



(4) 



(see Eqn. 2). Hereafter we take the characteristic size for th e high-ionization lines ( Hell 
Lyman complex) to be half this value, namely -Rblr — -Rblr/2 (jPoutanen fc Sternll2010l ). and 



fix ^BLR = 0.1. The opacity for the few-GeV photons is then Tyj{5 GeV) = r-j x (a^^/a-r), 
where a^^/a^ is the photon-photon pair production cross section in the units of the Thomson 
cross section, and 



5 GeV 



Hell Ly 



-2 (5) 

An -Rblr c Eueii 

for all the parameters as discussed above and the characteristic photon energy of the He II 
Lyman emission -Ehgii — 50 eV. Even though implying formally r^^ < 1 for -E ~ 5 GeV 7-ray 
photon^, this condition is exactly what is needed to produce a GeV break observed in the 
spectrum of 4C +21.35, because in the analyzed model one expects AF ~ r.|'^'^^/4 ~ 0.5. In 
other words, the fact that 4C +21.35 is characterized by a lower accretion disk luminosity 
than 3C 454.3 is consistent — in the framework of the discussed model — with the fact that 



Note that, assuming an isotropic distribution o f soft photons, the maximum of the pair production cross 



section cr^-y (Tt/3 is reached for Eq ~ im^c (jZdziarski 



198S 
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the spectral break in the former source is smaller than in the latter one (AF ~ 0.5 versus 



1/2 

1.2; note here that Lj^ja oc L^^sk and -Rblr cx: -Z^disk; ^^^^ tt oc L 



1/2 

disk; 



We conclude that as long as the Fermz-L AT o bserv ations of 4C +21.35 are considered 
solely, the model proposed by iPoutanen fc SternI (120101 ) can nicely account for the break 
observed in the spectrum of this blazar. An additional observational constraint is provided 
bv the detection of 4C +21.35 at 0.07 - 0.4 TeV energies by the MAGIC telescope jMariotti 
2010l : lAleksic et al.l [submitted to ApJLl ) , which coincided with the second GeV flare of the 
source. This coincidence suggests co-spatiality of the GeV and TeV emitting regions. Yet, if 
the GeV break is indeed due to the pair production on the high-ionization lines emitted by 
the BLR clouds, the emission zone has to be located at distances r ^ -Rblr — 10^'' cm from 
the central engine. At such distances, on the other hand, the photon-photon annihilation 
related to the near-infrared Aq — E^h/2mlc^ ~ 0.2 — 1 fim emission of the circumnuclear 
matter is likely to prevent the Ej ~ 0.07 — 0.4 TeV photons from escaping the circumnuclear 
environment, as argued below. 

First, let us note that the multi-body emission continuum of the standard accretion disk, 
peaking around ~ 10^^ Hz frequencies, extrapolates down to near-infrared (NIR) frequencies 
with the ~ 1 nm lumiriosity constituting at least 10% of the total disk luminosity (see, e.g., 
Tavecchio fc Ghisellinill2008l ). This, together with Ldisk estimated above, would then imply 
the NIR luminosity of the accreting matter in 4C +21.35 of the order of > 10^^ erg s~^, which 
is in good agreement with the valu e 9^,^ — 3.6 x lO^^ergs"^ following from the 2.2 /xm 
flux F2.2^m — 3.873 mJy provided by iFan et al.l (120061 ) for the object. Taking into account a 
possible contribution of the hot dust (located at some further distance, see below) and of the 
synchrotron jet emission to the observed NIR flux, we consider Lnir — lO^^ergs"^ to be a 
safe guess regarding the disk luminosity in 4C +21.35 indeed. Yet this emission is produced 
in the innermost parts of the accretion flow, with the characteristic size much smaller than 
-Rblr- This implies that the NIR photons emitted directly by the disk illuminate the jet 
region located around -Rblr exactly from behind, what reduces significantly the efficiency of 
the pair product ion process du e to the very kinematics of the photon-photon interaction (see 
the discussion in Reimer 2007 1^. Hence, the opacity for the < TeV radiation produced within 
the BLR is dominated not by the direct NIR emission of the accretion disk, but instead by 
the soft infrared photons emitted by the BLR itself, and also by the circumnuclear dust 
distributed beyond the BLR clouds. 



"'it may be easily demonstrated that for the jet regions placed much closer to the active center than -Rblr, 
this efficiency is not going to be drastically reduced. Therefore, a < TeV emission zone in 4C +21.35 located 
at r ^ -Rblr niay be excluded. 
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Regarding the former component , we re fer to the detailed calculations presented in 



Reimerl (120071 ) and iTavecchio fc MazinI (120091 ). which are based on the working assumption 



about a spherical distribution of BLR clouds around the active center. These calculations 
show in particular that the corresponding optical depth for the < 1 TeV photons is of the 
order a few, assuming the standard value ^blr ~ 0.1. This already suggests that the sub- 
TeV emission of 4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not 
likely to be produced inside the region of the highest ionization of BLR, but instead at 
r > -Rblr- Even more problematic may be however the latter component, which is related to 
the circumnuclear hot dust region (HDR) forming a torus-like structure around the accreting 
black hole. The characteristic scale of this structure, -Rhdr? can be estimated by applying 
the established relation between the disk lu minosity Lnisk and i?HDR (hereafter we assume 



the hot dust temperature Thdr — 10 K; see iSikora et al 



20091 . and references therein). For 



the particular parameters of 4C +21.35 one obtains in this way, roughly, 

^HDR ^ 4 ( -P^] PC ^ 10^9 cm . (6) 



1046 ergs-1 

The corresponding optical depth for the TeV photons produced within r < -Rhdr is then 

r^^(lTeV) ^ — — — 100 — — , (7) 

i/lTT itHDR C -&HDR V ^.i / 

where -Ehdr — 0.3 eV, and ^hdr is the fraction of the disk radiation reprocessed in the inner 
regions of the dusty torus, expected again to be at the level of 10%. 

The obtained value of t^7(1 TeV) is confusingly high. It may suggest that — z/ the 7-ray 
continuum of 4C +21.35 in the source rest frame extends up to ^ 1 TeV energies, while at 
the same time the emission zone is located at distances r < -Rhdr — the ^hdr parameter is 
significantly lower than expected (and, in fact, than inferred for the other analogous sources). 
Alternatively, the distribution of the hot dust around the nucleus of the discussed object 
may be very different from the simple geometry anticipated here (and in the literature). 
Nevertheless, even if both effects are relevant, the basic estimates presented above imply 
that one really should expect non- negligible, or even significant pair-production signatures 
around TeV photon energies in the intrinsic spectrum of 4C +21.35. Lack of such signatures 
would therefore place the blazar emission zone in the discussed objec t at relativel y large 



distances from the active center, r > -Rhdr, as recently advocated by ISikora et al.l (|2009[ ) 
for luminous blazars in general. On the other hand, the fact that the MAGIC spectrum of 
4C +21.35 extends up to 0.4 TeV energies only, may be providing a crucial piece of evidence 
in this context, since the emission of the hot dust is expected to fall sharply at wavelengths 
A < 1 /xm, thus decreasing the optical depth for < 0.4 TeV photons. 
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4. Conclusion 

In this paper we report on tlie two-year-long Fermi-LAT observation of the peculiar 
blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state from at least 
August 2008 until September 2009, and then entered a period of enhanced activity consisting 
of a gradual flux increase with some moderate flares superimposed. In April and June 
2010, 4C +21.35 underwent a very strong GeV outburst composed of several major flares 
characterized by the rise and decay timescales of the order of a day, and by the daily 7-ray 
fluxes nearly comparable to that of the Vela pulsar, the brightest persistent source in the 
7-ray sky. During the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law 
form with spectral breaks located near 1 — 3 GeV photon energies, similar to the flaring 
spectra of another luminous blazar, 3C 454.3. We demonstrate that, at least during the 
major flares, the jet in 4C +21.35 carries a total power comparable to the total accretion 
luminosity available to feed the outflow. As a result, the two major flares detected by 
the Fermi LAT should be considered as being saturated at the maximum level, and we 
should not expect to detect any flares from 4C +21.35 more luminous than -F>iooMeV ~ 
10~^phcm~^ s~^. We also discuss the origin of the break observed in the flaring spectra of 
4C +21.35. We show that, in principle, the model involving annihilation of the GeV photons 
on the He 1 1 Lyman recomb ination continuum and line emission of BLR clouds (as proposed 



by iPoutanen fc Stern! I2OIOI ) may account for such breaks. However, we also discuss the 
additional constraint provided by the surprising detection of 4C +21.35 at < ITeV energies 
by the MAGIC telescope, which coincided with the second GeV flare of the source. We argue 
that this sub-TeV emission (as well as the GeV emission of the source, if co-spatial), is not 
likely to be produced inside the region of highest ionization of BLR (~ 10^^ cm), because of 
the expected large opacity for 7-rays related to photon-photon annihilation on the infrared 
photon field provided by the accretion disk, BLR, and nuclear dust. Instead, it seems to 
originate further away from the active center, most likely as far as the characteristic scale of 
the hot dusty torus surrounding the 4C +21.35 nucleus (~ 10^^ cm). 

After completing this work, we learned that our basic estimates regarding the TeV 
opacity are supp orted by the most recent analysis of the Spitzer Space Telescope data by 



Malmrose et al.l (2011), who found a prominent infrared excess in the observed spectrum of 



4C +21.35, and showed that the excess is well modeled by the ~ 10^ K dust emiss ion. The 



extremely high isotropic luminosity of the hot dust estimated by iMalmrose et al.l indicates 
that the fraction of the disk radiation reprocessed in the inner regions of the dusty torus 
in 4C +21.35 is as large as ^hdr ~ 1- The exact spectral shape and the luminosity of the 
infrared continuum of the source available in a near future, as well as of the detail of the TeV 
spectrum soon reported by the MAGIC Collaboration (Aleksic et al., 2011, to be submitted), 
will allow us to re-examine the conclusions and estimates presented in our paper. 



- 16 - 



The Fermi-LAT Collaboration acknowledges generous ongoing support from a number 
of agencies and institutes that have supported both the development and the operation of the 
LAT as well as scientific data analysis. These include the National Aeronautics and Space 
Administration and the Department of Energy in the United States, the Commissariat a 
I'Energie Atomique and the Centre National de la Recherche Scientifique / Institut National 
de Physique Nucleaire et de Physique des Particules in France, the Agenzia Spaziale Italiana 
and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization 
(KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallen- 
berg Foundation, the Swedish Research Council and the Swedish National Space Board in 
Sweden. 

Additional support for science analysis during the operations phase is gratefully acknowl- 
edged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes 
Spatiales in France. 

Y. T. T. is supported by JSPS research fellowships for Young Scientists. L. S. is grateful 
for the support from the Polish MNiSW through the grant N-N203-380336. Y. T. T. and 
L. S. acknowledge Marek Sikora for his useful comments and remarks. The authors also 
thank A. Marscher for providing the results of the Spitzer observations of 4C -1-21.35. 

REFERENCES 

Abdo, A. A., et al. 2010a, ApJS, 188, 405 

— . 2010b, ApJ, 715, 429 

— . 2010c, ApJ, 723, 1082 

— . 2010d, ApJ, 716, 835 

— . 2009a, ApJ, 699, 817 

— . 2010e, ApJ, 710, 1271 

— . 2009b, ApJS, 183, 46 

Ackermann, M., et al. 2010, ApJ, 721, 1383 

Albert, J., et al. 2008, Science, 320, 1752 

Aleksic, J., et al. submitted to ApJL 



-17- 



Atwood, W. B., et al. 2009, ApJ, 697, 1071 

Carrasco, L., Carramifiana, A., Recillas, E., Porras, A., & Mayya, D. Y. 2010, The As- 
tronomer's Telegram, 2626, 1 

Celotti, A., & Ghisellini, G. 2008, MNRAS, 385, 283 

Celotti, A., Padovani, P., & Ghisellini, G. 1997, MNRAS, 286, 415 

Ciprini, S. 2009, The Astronomer's Telegram, 2349, 1 

Dominici, T., Pereyra, A., Abraham, Z., Almeida, U. B. D., Boisson, C, Chadwick, P., 
Magalhaes, A. M., & Ward, M. 2010, The Astronomer's Telegram, 2693, 1 

Donato, D. 2010, The Astronomer's Telegram, 2584, 1 

Fan, Z., Cao, X., & Gu, M. 2006, ApJ, 646, 8 

Finke, J. D., & Dermer, C. D. 2010, ApJ, 714, L303 

Ghiselhni, G., & Tavecchio, F. 2008, MNRAS, 387, 1669 

— . 2009, MNRAS, 397, 985 

Hartman, R. C, et al. 1999, ApJS, 123, 79 

Homan, D. C, Ojha, R., Wardle, J. F. C, Roberts, D. H., AUer, M. F., AUer, H. D., & 
Hughes, P. A. 2001, ApJ, 549, 840 

Hooimeyer, J. R. A., Schilizzi, R. T., Miley, G. K., & Barthel, P. D. 1992, A&A, 261, 5 

lafrate, G., Longo, F., & D'Ammando, F. 2010, The Astronomer's Telegram, 2687, 1 

Jorstad, S. G., Marscher, A. P., Mattox, J. R., Wehrle, A. E., Bloom, S. D., & Yurchenko, 
A. V. 2001, ApJS, 134, 181 

Kharb, P., & Shastri, P. 2004, A&A, 425, 825 

Kovalev, Y. Y., et al. 2005, AJ, 130, 2473 

Longo, F., Giroletti, M., & lafrate, G. 2009, The Astronomer's Telegram, 2021, 1 

Malmrose, M. P., Marscher, A. P., Jorstad, S. G., Nikutta, R., & Elitzur, M. 2011, submitted 
to ApJ 

Mariotti, M. 2010, The Astronomer's Telegram, 2684, 1 



- 18 - 



Mattox, J. R., et al. 1996, ApJ, 461, 396 
Moderski, R., & Sikora, M. 2005, Ap&SS, 297, 369 

Nandikotkur, G., Jahoda, K. M., Hartman, R. C, Mukherjee, R., Sreekumar, P., Bottcher, 
M., Sambruna, R. M., & Swank, J. H. 2007, ApJ, 657, 706 

Neronov, A., Semikoz, A. N. D., & Vovk, I. 2010a, The Astronomer's Telegram, 2617, 1 

Neronov, A., Semikoz, D. V., & Vovk, I. 2010b, ArXiv e-prints 

Nesci, R., & Montagni, F. 2010, The Astronomer's Telegram, 2708, 1 

Osterbrock, D. E., & Pogge, R. W. 1987, ApJ, 323, 108 

Poutanen, J., & Stern, B. 2010, ApJ, 717, L118 

Reimer, A. 2007, ApJ, 665, 1023 

Reynolds, S. P. 2009, ApJ, 703, 662 

Saikia, D. J., Wiita, P. J., & Muxlow, T. W. B. 1993, AJ, 105, 1658 
Sikora, M., & Madejski, G. 2000, ApJ, 534, 109 

Sikora, M., Madejski, G., Moderski, R., & Poutanen, J. 1997, ApJ, 484, 108 
Sikora, M., Stawarz, L., & Lasota, J. 2007, ApJ, 658, 815 

Sikora, M., Stawarz, L., Moderski, R., Nalewajko, K., & Madejski, G. M. 2009, ApJ, 704, 38 
Striani, E., et al. 2010, The Astronomer's Telegram, 2686, 1 
Tavecchio, F., & GhiseUini, G. 2008, MNRAS, 386, 945 

Tavecchio, F., GhiseUini, G., Bonnoh, G., & Ghirlanda, G. 2010, MNRAS, 405, L94 
Tavecchio, F., & Mazin, D. 2009, MNRAS, 392, L40 
Verrecchia, F., et al. 2010, The Astronomer's Telegram, 2698, 1 
— . 2009, The Astronomer's Telegram, 2348, 1 

Wagner, S. J., Behera, B., & HESS collaboration. 2010, in Bulletin of the American Astro- 
nomical Society, Vol. 41, Bulletin of the American Astronomical Society, 660 — h 

Wang, J., Luo, B., & Ho, L. C. 2004, ApJ, 615, L9 



- 19 - 



Zdziarski, A. A. 1988, ApJ, 335, 786 



This preprint was prepared with the A AS macros v5.2. 



-20- 



State 



Start date (MJD) 



End date (MJD) 



Quiescent 

Intermediate+ Active 
Intermediate 
First flare 
Inter-flare 
Second flare 



2008 August 4 (54682.66) 2009 September 12 (55086.00) 



2009 September 12 (55086.00) 
2009 September 12 (55086.00) 

2010 April 23 (55309.00) 
2010 May 2 (55318.00) 

2010 June 16 (55363.00) 



2010 June 24 (55371.00) 
2010 April 23 (55309.00) 
2010 May 2 (55318.00) 
2010 June 16 (55363.00) 
2010 June 24 (55371.00) 



Table 1: Time intervals defined in this paper; see Figure 1 and 2 (a). 



State 




The 




AL 








[GeV] 




Intermediate-I- Active 


2.18±0.02 


2.64±0.06 


2 4+°-2 


-27.2 


Intermediate 


2.30±0.05 


2.69±0.40 




-6.6 


First flare 


1.80±0.06 


2.40±0.07 


-| 1+0.3 


-17.9 


Inter-flare 


2.24±0.03 


2.81±0.14 


9 7+0.6 
^- '-0.6 


-8.4 


Second flare 


2.00±0.05 


2.44±0.10 


-•-■'-0.4 


-6.3 



Table 2: Spectral fits using broken power-law model with low- and high-energy spectral 
indices Fle and Fhe, respectively, and the break photon energy E^,-. In the last column, AL 
is the difference of the logarithm of the likehhood of the fit with respect to a single power-law 
fit. 



State 


a 


P 


AL 


Intermediate+Active 


2.11±0.03 


0.06±0.01 


-19.6 


Intermediate 


2.26±0.05 


0.04±0.02 


-3.1 


First flare 


1.70±0.09 


0.13±0.02 


-17.9 


Inter-flare 


2.14±0.06 


0.07±0.02 


-7.0 


Second flare 


1.94±0.08 


0.07±0.02 


-4.8 



Table 3: Alternative spectral fits using log-parabola function dN/ dE oc {E / Ecr)~"~ " i 

where E^r is fixed to 300 MeV. In the last column, AL is the difference of the logarithm of 
the likelihood of the fit with respect to a single power-law fit. 
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Fig. 1. — Weekly (7-day-bin) 7-ray light curve of 4C +21.35 (photon energies E > 100 MeV) 
from 2008 August to 2010 August obtained by the Fermi LAT. The 95% upper limits are 
shown by red triangles. Three main epochs marking different activity levels of the source 
are defined (see also Table 1 ) . 
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Fig. 2. — (a) Daily 7-ray light curve {E > 100 MeV) during an active phase of 4C +21.35 
as observed by the Fermi LAT. Black arrows indicate the epochs when two very high en 
ergy p hotons with energies above 100 GeV were detected by the Fermi LAT ( iNeronov et al. 
2010aJ). A red arrow depic ts the date MJ D 55364.9 when MAGIC detected very high energy 
emission from 4C +21.35 ( lMariottill2010l ). (b) Daily variation of the photon index derived 
from a single power-law fitting to the LAT data, (c) Relation between daily 7-ray fluxes 
and photon indices measured during the first flare. The numbers indicate the days from the 
onset of the flare, (d) The same as (c) but for the second flare. 
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Fig. 3. — The spectral energy distribution of 4C +21.35 within the LAT energy range during 
each activity state of the source defined in Figure 1 and 2 (a) (see also Table 1). Solid lines 
show the best-fit broken power-law models. Arrows indicate 95% upper limits for the binned 
spectra. 
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Fig. 4. — The three different spectral energy distributions of 4C +21.35 corresponding to 
different periods of the active state of 4C +21.35 (red, blue, and green filled circles), compared 
with the flaring spectrum of 3C 454.3 ('the bright D ecember 2009 flare spec trum' of 3C 454.3 
corresponding to the epoch MJD 55166-55173; see lAckermann et al.ll2010l gray diamonds). 
Solid lines show the best-fit broken power-law models. Arrows indicate 95% upper limits. 



